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An overshoot solar dynamo with a strong return meridional flow 
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The meridional circulation plays an essential role in determining the basic mechanism of the dynamo action in case the of a low eddy 
diffusivity. Flux-transport dynamos with strong return flow and a deep stagnation point are discussed in the case of a positive «-effect 
located in the overshoot layer and a rotation law consistent with helioseismology. By means of a linear dynamo model, it will be shown 
that the migration of the toroidal belts at lower latitudes and the periods of the activity cycles are consistent with the observations. 
Moreover, at variance with previous investigations, the typical critical dynamo numbers of dipolar solutions are significantly smaller that 
those of quadrupolar solutions even in the regime of strong flow. 
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1 Introduction 

The flux-transport dynamo is a promising mechanism to explain several properties of the solar activity 
cycle. Although some aspects still need to be clarified, the basic observation behind this process is that in 
the presence of a low eddy diffusivity r]t the magnetic Reynolds number LU/r]t becomes very large and 
the dynamics of the mean- field flow U is thus an essential ingredient of the dynamo process. In this regime 
the advection produced by the meridional circulation dominates the diffusion of the magnetic field which 
is then "transported" by the meridional circulation. Clearly the word "transport" must be used with some 
care because the field is not completely "frozen" into the plasma: rather, the propagation of the dynamo 
wave is significantly distorted so that the well-known Parker- Yoshimura law is not necessarily satisfied in 
this situation. 

Local helioseismology and doppler speed measurements substantially agree on a detection of an av- 
erage surface fiow of about 15 ms~^ around 30°, with a peak of about 20 — 25 ms~^. The return 
flow, located near the base of the convection zone, is much more difficult to detect and only theoret- 
ical estimates are at our disposal. In fact previous models of flux-transport dy namo have cons i dered 



a return flow of about one order of magnitude sinaller than the surface flow (iBonanno et al.\ |2002| . 
Dikpati and Charbonneaull 19991, Ch atterie e et adl2004 . Guerrero and de Gouveia Dal Pin 0I I2OO8I I although 
in recent studies ( Kiiker and Riidig er .. ,20081 . iKiiker et a/.ll20li r it has been argued that the strength of this 
flow is instead of the same order of the poleward surface flow. 

The other important ingredients of the model are the a-effect and the differential rotation. While the 
latter can be deterr nined by helioseismology, for t he former it will be assumed that the most suitable loca - 
tion for this effect ( Parker 19551 . Steenbeck et al\ \l966) is just ben eath the convection zoii e (|Parkerlll993l ;) 
where a strong radial shear is produced in the so called tachocline ( Spiegel and Zahn 19921 ). The source of 
the turbulent helicity producing the a-effect can be attributed to variou s mechanisms. The most promising 
one is the tachocline instability proposed by Dikpati and Gilman ( 200ll ). although in recent investiga tions 
another appealing possibility is provided by a current-helicity generated a- effect (iGellert et ai.ll200lh due 
to kink and quasi-interchange instabilities in stably stratified plasmas ( Bonanno and Urpin " l201ll ). The 
aim of this work is to study models of fiux-transport dynamo with a strong meridional fiow, where the 
stagnation point is self-consistently determined fr om the balance of the angular momentum in a turbulent 
plasma, following the self-consistent approach of Durney (j2000l ). It will be shown that dipolar solutions 
are strongly favored and that, in the advection dominated regime, the period is essentially determined 
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by the strength of the retu rn flow, at variance with the conclusion obtained by a recent investigation 
( Pipin and Kosovichevl [2OI1I ) . 



2 Basic Equations 

As is weh known, the magnetic induction equation reads 

Vx{U xB + aB) - V X {rjtV x B) 



dB 

'dt 



(1) 



where rjt is the turbulent diffusivity. Axisymmetry implies that relative to spherical coordinates the mag- 
netic field B and the mean fiow field V are given by 



B = B^r, e, t)e^ + V x [A{r, 6, t)e^ 



V = u{r,9) + r sin 0Q{r, 6)e^, 



(2) 



where B(f,(r, 6, 4>)e^ and V x [A{r, 6, t)e(f,] are the toroidal and poloidal components of the magnetic field re- 
spectively. Moreover, the meridional circulation u(r, 6) and differential rotation I7(r, 9) are the poloidal and 
toroidal components of the global velocity flow fleld V. In particular the poloidal and toroidal components 
of ^ respectively determine 
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where s = rsin^. The a-effect is always antisymmetric with respect to the equator, so that we write 



a = - ao cos 6* 



1 + erf 



ai 



d 



1 — erf 



X - a2' 
d . 



(4) 



where ao is the amplitude of the a-effect, x = v/Rq is the fractional radius, ai, 02 and d deflne the location 
and the thickness of the turbulent layer. In our investigation, we use the values ai = 0.68, 02 = 0.72 and 
d = 0.025. It is reasonable to imagine that below the tachocline the turbulent diffusivity decreases by a 
few orders of magnitude from the turbulent value attained in the bulk of the convection zone. We can 
conveniently represent this transition with the following functional form 



I, 

ric + ^{rit 



1 + erf 



(5) 



where r]t is the eddy diffusivity, r]c the magnetic diffusivity beneath the convection zone and d^i represents 
the width of this transition. In particular, we use the values rjt/rjc = 10^, d = 0.02 and = 0.71. 

The components of the meridional circulation can be represented with the help of a stream function 
lF(r, e) = — sin^ cos e ip{r) so that 



1 d^ 1-3 cos2 e 



r'^psinO de 



pr^ 



ip{r), 



U0 



1 dW cos sin dtp{r) 



rpsm.e dr 



pr 



dr 



(6) 



with the consequence that the condition 'V-{pu) = is automatically fulfllled. A positive ip describes a cell 
circulating clockwise in the northern hemisphere, i.e. the flow is polewards at the bottom of the convection 
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Figure 1. The a-effect (solid line), turbulent diffusivity (dot-dashed line) and (minus) the function ipif) (dashed line), are depicted 
in the left panel. The streamlines of the flow are shown in the middle panel, dashed lines label the anticlockwise flow. The meridional 
circulation in units of the maximum surface value at a latitude of 45° is instead depicted in the right panel. 



zone and equatorwards at the surface. For a negative the flow is, as is observed, polewards at the surface. 
In order to keep the flow inside the convection zone, the function -0 must be zero at the surface and at the 
bottom of the convection zone. The helioseismic profile for the differential rotation is taken so that 



f2(r,0) = f2e + ^[l+erf(^^) - r?,) , (7) 



where ]7c/2vr = 432.8 nHz is the uniform angular velocity of the radiative core, f2s{6) = Q^q + ^2 cos^ 9 + 
04 cos^ 9 is the latitudinal differential rotation at the surface. In particular ]?eq/2'7r=460.7 nHz is the angular 
velocity at the equator, a2/2vr=— 62.9 nHz and a4/27r=— 67.13 nHz, erf(x) is the usual error function. In 
this calculation the angular velocity is normalized in terms of equator differential rotation fieq-, = 0.71 
and dc = 0.025. As usual, the dynamo equations can be made dimensionless by introducing the dynamo 
numbers 

— ) — ) ^uj — ! — ) \P) 

m vt m % 

where R is the stellar radius, lo is the frequency of the dynamo wave and U = ug{r = R,9 = 45°). The 
meridional circulation is largely unknown, although helioseismology can provide an upper limit to the 
strength of the poleward flow as discussed in the introduction. A strategy to constraint several properties 
of the meridional circulation is to assume the differential rotation profile i7(r, 9) as a given ingredient, and 
deduce an approximat ion for the func tion from the angular momentum conservation along the azimuthal 
direction as shown by Durnev ( 2000l ). An approximate expression for ■(/; is thus 



^~ ^ [\urue)d9 (9) 

in which r the Coriolis number. In particular, for the standard, isotropic mixing-length theory ([9]) becomes 
(see Durney 2000 for details) 

il) ^ -f>prT{ul) . (10) 

In principle it would be possibile to explicitly compute ij^ and tp' using the relation (llOp knowing the 
convective velocities of the underlying stellar model. In practice this would be problematic, because the 
convective fluxes and their radial derivatives computed from standard MLT are discontinuous at the base 
of the convective zone. In a more realistic situation the presence of an overshoot layer implies that (u^) — )• 
smoothly so that uq is continuous at the inner boundary. Nevertheless one can use the representation (jlOp 
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to determine the stagnation point where '0' = 0, which turns out to be around r = 0.8 solar radii in a 
standard solar model. 

This value is only slightly greater than the value obtained by iKiiker et al. where the stagnation 



point is at r = 0.77 (private communication). On the the other hand those authors uses a stress- free 
boundary condition for the meridional circulation that implies a non-zero value of the flow at the inner 
boundary. The problem is that in the overshoot region the eddy diffusivity drops of several orders of 
magnitude and the Reynolds number becomes very large at the inner bottom, most probably producing 
a boundary layer instability. I argue that stress free boundary conditions for the meridional circulation in 
the context of flux-transport dynamo are probably unphysical and one should try to consider a situation 
where both Ur and ug vanishes at the inner boundary. 
An explicit form of the function ^ which incorporates the following features reads 



= c 



exp 



{X - Xfc) 
C72 



(11) 



where C is a normalization factor, Xf, = 0.67 defines the penetration of the flow, a = 0.025 measures how 
fast (n^) decays to zero in the overshoot layer and the location of the stagnation point. The density profile 
is taken to be 

/I \ 

P = Po(--xo) (12) 

in which m is an index representing the the stratification of the underlying solar model, its value in the 
region of interest is approximately 2, and xq = 0.9. The radial profile of the a-effect, turbulent diffusivity, 
stream function and meridional circulation used in the calculation is depicted in figure ([I]). 

This linear dynamo problem is solved with a finite-difference scheme for the radial dependence and a 
polynomial expansion for the angular dependence. In particular, the following expansions for the field are 
used: 

A{r, e) = e^* a„(r) P^(cos 9), (13a) 

n 

B{r, 9) = e^* J] &„(r) P^(cos 9), (13b) 

m 

where A is the (complex) eigenvalue so that Im{A} = u, the frequency of the dynamo wave, n = 1,3,5,... 
and m = 2,4,6,... for antisymmetric modes, and n m for symmetric modes. Vacuum boundary 
conditions at the surface are then translated into 

^ + (n + 1) a, = 6^ = 0. (14) 

In the interior at x = Xi = 0.6 we have instead the set 

x-^ + hm = an = ^, (15) 

which imply perfect conductor boundary conditions. 

On substituting (I13h .b) into ([3K,b) one obtains an infinite set of ordinary differential equations that 
can be conveniently truncated in n when the desired accuracy is achieved. The system is in fact solved 
by means of a second order accuracy finite difference scheme and the basic computational task is thus to 
numerically compute eigenvalues and eigenvectors of a block-band diagonal real matrix of dimension Mxn, 
M being the number of mesh points and n the number of harmonics, M{a)v = Xv and v is in general 
a complex eigenvector. This algorithm is embedded in a bisection procedure in order to determine the 
critical Co-value needed to find a purely oscillatory solution, for which Re{A} = 0. For actual calculation 
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Figure 2. Butterfly diagram for a overshoot dynamo solution with no meridional circulation and two different angular dependence of 
the a-effect. In the left panel the a-effect is proportional to cos0 and the solution has Ca = 19.34, Cn = 3000, and C'u = 33.22. The 
solid and dashed lines represent the radial field at the surface (solid for negative Br), and blue is for negative toroidal field. The right 
panel represents instead a model with an a-effect proportional to cos9 sin^ 9 with Ca = 15.34 and Cui = 23.02. (colour online) 

the resolution of 60 x 60 has been used because it was checked that further terms in the polynomial 
expansion and in the mesh points did no t lead to any significant change in the solution. The code has been 
extensively tested in lJouve et al. (l2008l ^. 



3 Results 

In this section the effect of a strong return flow with a deep stagnation point will be discussed. It is 
instructive to first consider the case of vanishing meridional circulation, where the propagation of the 
dynamo wave is determined by the Parker- Yoshimura law. In particular if the a-effect has the simple cos 9 
angular dependence the generation of the poloidal field from the toroidal is more pronounced at higher 
latitudes and, as dfi/dr < there, the dynamo wave propagates towards the equator until it reaches the 
latitude where dfi/dr « and the migration is halted. The solution is characterized by the following 
dynamo numbers, Cq, = 19.34, Cq = 3000, and frequency = 33.22, which implies a period of about 6.2 
years for the Sun given a turbulent diffusivity of rjt = 4.66 lO^'^cm^s"^. Models with a negative a-term at 
the bottom of the convection zone produce only stationary solutions. The butterfly diagram is depicted in 
the left panel of flgure (l2|), where the radius is taken at maximum B^. It is interesting to see what happens 
for an angular dependence of the type cos 6 siv? 6 for the a-term since the field regeneration due to the 
dynamo action occurs at lower latitude in this case. The solution for a positive a-term is depicted in the 
right panel of figure ([2]). In this case the dynamo action occurs around latitudes at which dQ/dr ~ and 
two distinct branches, one equatorward caused by the dQ/dr < at high latitudes and the other poleward, 
at lower latitude, caused by dQ/dr < are present. The poleward migration of the butterfly diagram at 
low latitude and the phase relation are clearly not consistent with the observations. In the case of a negative 
a-term the only possible solution is a stationary one, as before. The conclusion of the above calculations 
is that it is difficult to imagine that a simple overshoot dynamo can successfu lly reproduce the observed 



featu res of the solar cycle, difficulties also noticed in earlier investigations (iBrandenburg and Riidiger 
1995h 



The situation radically changes in the flux-dominated regime, ^ 1 as one can see in figure ([3]) as in 
this case the meridional circulation transports the flux at lower latitudes. This solution is characterized 
by a flow U = 27 ms^^ with rjt = 4.66 lO^^cm^s^^ so that C„ = 400. The relevant question is to 
discuss the parity of the solutions as a function of the strength of the meridional circulation since a strong 
equatorward flow at the bottom of the convection zone can easily favor a quadrupolar parity. In fact in 
this case the field is by definition non-zero at the equator and the presence of the meridional circulation 
can help the occurrence of solutions with this symmetry . Usu ally only a small window in the parameter 
space is consistent with the observations ( Bonanno et al. 20061 ). 



However, as it is shown in figure ([5]) the critical Ca value is significantly smaller for dipolar solutions 
(solid line) than for quadrupolar solutions (dashed line), at variance with previous investigations where 
the difference in the critical dynamo numbers between dipolar and quadrupolar solution was not found 






Figure 3. The temporal evolution of the magnetic field for a solution with Cu = 400, Ca = 19.34, Cq = 30000, and frequency Cuj = 102 
which imply a period of about 20 yrs and a flow of 27ms~^. The left part are the isocontours line of the toroidal field, with solid line 
for negative and dashed line for positive value of the field. The right part represents of the streamlines of the poloidal field given by 
contours of ArsinS. Solid line are for negative values of A. 



to be very significative ( Bonanno et al. 20021 ) . It is important to stress that the evolution of the critical 
dynamo number Ca and the period is not a monotonic function of the meridional circulation, in general. 
As it is clear from figure (j5]) there are essentially two different scaling regimes for weak and for strong flow. 
In the first case the meridional circulation does not significantly affect the behavior of the dynamo wave, 
while at large value of the fiow, the period decreases as the flow increases. The two scaling regimes are 
separated by a region of stationary solution around 7ms~^ for an eddy difFusivity of rit = 4.6 6 lO^^cm^s"^. 
This result is in agreement with the flndings of previous studies (jBonanno et a/Jl2002l . [2OO6I ) and it is not 
changed by th e presence of a strong flow wi th a deep stagnation point. I conclude that the difference with 
the flndings of Pipin and Kosovichev ( 20 111 ) are mainly due to the inclusion (in that work) of the f2 x J 
contribution in the turbulent electromotive force (see their figure 5 where a saturation of the period seems 
to occur instead). 
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Figure 4. Butterfly diagram and phase relation for the solution presented in figure l(3j- (colour online) 
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Figure 5. The critical dynamo numbers (left panel) and the period (right panel) of the solutions as a function of the meridional circulation 
for dipolar solutions (solid line) and quadrupolar solution (dashed lines) for an eddy difEusivity of rji = 4.66 lO^^cm^s"^ 

4 Conclusions 



The results presented in this investigation suggest that a standard a^O mean field dynamo with a deep 
meridional circulation and an a-efFect located in the overshoot layer reproduces the migration of the 
toroidal belts at lower latitude in dynamo models with heliosesimically derived rotation law. At variance 
with previous studies (with the exception of the recent work by Pipin & Kosovichev 2011), the return flow 
is of the same order of the surface flow, and the solutions are clearly dominated by models with dipolar 
parity. 

In this investigation the super critical dyna i no act ion has not been considered. This point was investigated 
(for a different flow pattern) in ljouve et al\ 1 20081 ) where it was shown that values of Ca about ten times 
supercritical does not significantly change the nature of the solution, even in the flux-transport regime. 
On the other hand, the difference between dipolar and quadrupolar solutions in the fiux-dominated regime 
are rather significant, and one could argue that at least for not too large supercritical dynamo numbers 
the conclusions of this investigations will still be valid. 

The idea of determining the sta gnation point o f the flow by the maximum of the convective fluxes in the 
convection zone, as sugg ested bv IPurnevI (|200d l. leads to successful models of the meridional circulation 
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in flux-dominated dynamo. It would be important to extend the findings of this work to the stellar case, 
at least for slowly rotating solar-like stars. 
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